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Using 7.3 fb _1 of pp collisions collected by the DO detector at the Fermilab Tevatron, we measure 
the distribution of the variable tf>^, which probes the same physical effects as the Z/y* boson 
transverse momentum, but is less susceptible to the effects of experimental resolution and efficiency. 
A QCD prediction is found to describe the general features of the 0* distribution, but is unable to 
describe its detailed shape or dependence on boson rapidity. A prediction that includes a broadening 
of transverse momentum for small values of the parton momentum fraction is strongly disfavored. 

PACS numbers: 12.38.Qk, 13.85.Qk, 14.70.Hp 



Z/j* bosons are produced at hadron colliders via 
quark-antiquark annihilation. Their decays to e + e _ and 
can be detected with little background and the 
phenomenology is simplified by the absence of color flow 
between the initial and final states, thus providing an 
excellent testing ground for QCD predictions. Resum- 
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mation techniques |l| allow calculations of the distribu- 
tion of Z/j* boson transverse momentum, p^, within 
the framework of perturbative QCD, even at relatively 
low p" (e.g., < 30 GeV). However, additional non- 
perturbative form factors must be determined in global 
fits to experimental data [2j . An increase of these form 
factors for x < 10 -2 , where x is the parton momentum 
fraction, was suggested [3j to improve the description of 
hadron production observed in deep inelastic electron- 
proton scattering at HERA. Since vector boson produc- 
tion corresponds typically to parton x < 10~ 2 at the 
LHC, these modified form factors would lead to a broad- 
ening of the expected vector boson transverse momentum 
distributions [J] . This "small- a; broadening" would influ- 
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ence the measurement of the W boson mass as well as 
searches for Higgs bosons and physics beyond the stan- 
dard model at the LHC. It is important to study quanti- 
tatively such x-dependencies at the Tevatron, where they 
can be probed using the dependence of the p^ distribu- 
tion on boson rapidity 

In the region of low p^, the precision of the most 
recent measurements at the Tevatron 0, @] was domi- 
nated by uncertainties in correcting for experimental res- 
olution and efficiency. Furthermore, the choice of bin 
widths was restricted by experimental resolution rather 
than event statistics. The variable ax, which corresponds 
to the component of pff that is transverse to the dilep- 
ton thrust axis, t, has been proposed as an alternative 
analysing variable that allows us to study the issues dis- 
cussed above, but is less susceptible than the pj, to detec- 
tor effects Q. Figure [1] illustrates this and other relevant 
variables defined below. The ax distribution was subse- 
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FIG. 1: Illustration of the variables defined in the text and 
used to analyse the dilepton transverse momentum. 

quently calculated to ncxt-to-lcading-log (NLL) accuracy 
using resummation techniques [§]. Additional analysing 
variables with even better experimental resolution have 
recently been proposed and studied 10]. The optimal 
variable was found to be <j)*, which is defined as: 

4>* = tan(<£ acop /2)sin(0*), 

where acop is the acoplanarity angle, given by: acop = 
7r — A<fi u , and A(f> u is the difference in azimuthal an- 
gle, 4>, between the two lepton candidates. The vari- 
able 9* is a measure of the scattering angle of the lep- 
tons with respect to the proton beam direction in the 
rest frame of the dilepton system. It is defined 
by: cos(#*) = tanh \(wT — V + ) /2], where ry~ and i] + arc 
the pseudorapidities 5] of the negatively and positively 
charged lepton, respectively. 

The variable <fi* is highly correlated with the quan- 
tity ar/mee, where mu is the dilepton invariant mass. 
Since acO p and 6* depend exclusively on the directions 
of the two leptons, which are measured with a precision 
of a milliradian or better, 0* is experimentally very well 
measured compared to any quantities that rely on the 
momenta of the leptons. 



We present a measurement of the normalized <t>* dis- 
tribution, (1/er) x (da/dcf)*), in bins of \y\, using 7.3 fb _1 
of pp collisions collected by the DO detector at the Fer- 
milab Tevatron. The (f>* distributions are measured in 
both dielectron and dimuon events and are corrected for 
experimental resolution and efficiency. We correct back 
to the level of observable, generator-level particles; that 
is, we apply kinematic selection criteria at the particle 
level that match those applied in the selection of can- 
didate events in the data ll|. Particle level electrons 



are defined as the four- vector sum of any electrons and 
photons within a cone of ATI = ^{Ar/) 2 + (A0) 2 < 0.2 
around an electron, where A77 (A(f>) is the distance in 
r\ ((f>) from the particle level electron; this mimics the 
measurement of electron energy in the calorimeter. Par- 
ticle level muons are defined after QED final state radi- 
ation; this mimics the measurement of muon momentum 
in the tracking detector. The kinematic selection criteria 
are: electrons must satisfy px > 20 GeV and \rj\ < 1.1 
or 1.5 < \rj\ < 3; muons must satisfy px > 15 GeV and 
\rj\ < 2; mu must fall within the range 70-110 GeV. 

The corrected data are compared to predictions from 
the Monte Carlo (MC) program ResBos [12| with the 
above kinematic selection criteria applied at the particle 
level. ResBos generates Z/j* boson events with initial 
state QCD corrections to next-to-leading order (NLO) 
and NLL accuracy together with: a non-perturbative 
form factor, whose width is controlled primarily by the 
parameter g 2 (with default value [0.68jl ; 1 l G eV 2 )|2j; an 
additional next-to-NLO (NNLO) K-factor [H; CTEQ6.6 
NLO parton distribution functions (PDFs) [Hj]; and 
QED radiative corrections from PHOTOS [H| . The QCD 
factorization and renormalization scales are set event by 
event to the mass of the Z/j* boson propagator. 

The DO detector [lj| consists of: silicon microstrip 
and central fiber tracking detectors, located within a 2 T 
superconducting solenoid; a liquid-argon/uranium sam- 
pling calorimeter; and an outer muon system consisting 
of tracking and scintillation detectors located before and 
after 1.8 T toroids. Candidate dielectron events are re- 
quired to satisfy a trigger based on the identification 
of a single electron and to contain two clusters recon- 
structed in the calorimeter with a transverse and longi- 
tudinal shower profile consistent with that expected of 
an electron. The calorimeter is housed in three sepa- 
rate cryostats; this has the effect that electron identi- 
fication is degraded in the region 1.1 < |ry| < 1.5. Can- 
didate dimuon events are required to satisfy a trigger 
based on the identification of a single muon and to con- 
tain two muons reconstructed either in the outer muon 
system, or as an energy deposit consistent with the pas- 
sage of a minimum-ionizing particle in the calorimeter. 
In order to ensure an accurate measurement of the lep- 
ton directions at the point of production, the two lep- 
ton candidates are required to be matched to a pair of 
oppositely charged particle tracks reconstructed in the 
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central tracking detectors. Candidate leptons resulting 
from misidentified hadrons or produced by the decay of 
hadrons are suppressed by requiring that they be isolated 
from other particles in the event and, in the case of elec- 
trons with 1 77 1 < 1.1, by requiring the energy measured in 
the calorimeter and the momentum measured in the cen- 
tral tracking detectors to be consistent. Contamination 
from cosmic ray muons is strongly suppressed by a re- 
quirement that the muons originate from the pp collision 
point and by rejecting events in which the two muon can- 
didates are back to back in r\. In total, 455k dielectron 
events and 511k dimuon events are selected. 

The corrections to the observed <fi* distribution for 
experimental resolution and efficiency are evaluated us- 
ing Z/j* boson MC events that are generated with 
PYTHIA [ijj and passed through a GEANT-based [HI sim- 
ulation of the detector. These fully simulated MC events 
are re- weighted at the generator level in two dimensions 
(Pt and \y\) to match the predictions of ResBos. In ad- 
dition, adjustments are made to improve the accuracy of 
the following aspects of the detector simulation: electron 
energy and muon px scale and resolution; track </> and 
•q resolutions; trigger efficiencies; and relevant offline re- 
construction and selection efficiencies. Variations in the 
above adjustments to the underlying physics and the de- 
tector simulation are included in the assessment of the 
systematic uncertainties on the correction factors. 

The systematic uncertainties due to electron energy 
and muon pt scale and resolution are small, and arise 
only due to the kinematic requirements in the event se- 
lection. The measured 4>* distribution is, however, sus- 
ceptible to modulations in (f> of the lepton identification 
and trigger efficiencies, which result, e.g., from detector 
module boundaries in the calorimeter and muon systems. 
Particular care has been taken (a) in the choice of lepton 
identification criteria in order to minimize such modula- 
tions and (b) to ensure that such modulations are well 
simulated in the MC. For example, the requirements im- 
posed on shower profile are much looser than those usu- 
ally employed in electron identification within DO, be- 
cause tight requirements are particularly inefficient in the 
regions close to module boundaries in the calorimeter. 
Similarly, the inclusion of muon candidates identified in 
the calorimeter reduces the effect of gaps between mod- 
ules in the outer muon system. Accurate modelling of 
the angular resolution of the central tracking detectors 
is another crucial aspect of this analysis. The resolu- 
tion in </> and r\ is measured in the data using cosmic ray 
muons that traverse the detector, since these should pro- 
duce events containing two tracks that are exactly back 
to back except for the effect of detector resolution. 

Backgrounds from Z — ► t~t + , W — >• Iv (+jets), and 
WW — > Ivtv are simulated using pythia. Background 
from top quark pair events is simulated with ALPGEN [l!| , 
with pythia used for parton showering. Background 
from multijet events is estimated from data. The total 



fraction of background events is 0.26% for the dielectron 
channel, and 0.38% for the dimuon channel. 

Since the experimental resolution in <$>* is narrower 
than the chosen bin widths, the fractions of accepted 
events that fall within the same bin in <f>* at the par- 
ticle level and reconstructed detector level in the MC 
are high, having typical (lowest) values of around 98% 
(92%). Therefore, simple bin- by-bin corrections of the 
4>* distribution are sufficient. In almost all (/>* bins the 
total systematic uncertainty is substantially smaller than 
the statistical uncertainty. 

Figure [3] shows the corrected particle level (j>* distri- 
butions together with predictions from ResBos. Fig- 
ure [3] shows the ratio of the corrected </>* distributions 
to the ResBos predictions in both the dielectron and 
dimuon channels. The general shape of the distributions 
is broadly described by ResBos over the full range in 
<fi*. However, the small statistical uncertainties resulting 
from the large dilepton data sets, combined with the fine 
binning and small systematic uncertainties resulting from 
the use of <fi* as the analysing variable, reveal differences 
between the data and ResBos. Since the particle level 
definitions for electrons and muons to which the data 
are corrected are slightly different, Fig. [3] represents the 
most appropriate way to demonstrate the consistency of 
the dielectron and dimuon data. Given that the exper- 
imental acceptance corrections are very different in the 
two channels, this consistency represents a powerful cross 
check of the corrected distributions. 

The results of fits for the value of <?2, separately in each 
\y\ bin and channel, are shown in Tabic [I] It can be seen 
that the fitted values of g-2 show a monotonic decrease 
with increasing \y\ for both channels. That is, the width 
of the <f>* distribution becomes narrower with increasing 
\y\ faster in the data than is predicted by ResBos. This 
is the opposite of the behavior expected from the small- 
s' broadening hypothesis 0, Q • Figure [3] confirms that 
the prediction from ResBos with small- a; broadening is 
in poor agreement with data. It can also be seen that 
choosing the 172 value (0.66 GeV 2 ) that best describes 
the average behavior of the data over all \y\ bins and 
channels has little effect on the level of agreement with 
data. 



Channel 


\y\ < 1 


K\y\ <2 


|V|>2 


ee 


0.644 ± 0.013 


0.619 ± 0.017 


0.550 ± 0.048 




0.670 ± 0.012 


0.645 ± 0.019 





TABLE I: Value of g 2 (GeV 2 ) that best describes the cor- 
rected data in each \y\ bin and channel. 

A previous measurement 0) showed that, for central 
rapidities, ResBos underestimates the number of Z/7* 
bosons at high p!^ by about 10%. This is consistent with 
the deviations seen at high values of </>* in Fig. [3] (a). 

In summary, using 7.3 fb _1 of pp collisions collected 
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D0 7.3 fb" 1 
» ,li,u data 

ee data 
ResBos 

ResBos (small-x) 




FIG. 2: (color online) Corrected distributions of (1/cr) x (da/d(j>*) for dimuon events with (a) \y\ < 1 and (b) 1 < \y\ < 2; 
and dielectron events with (c) \y\ < 1, (d) 1 < \y\ < 2 and (e) \y\ > 2. The larger plots show the restricted range 
< (j>^ < 0.34 and the insets show the full range of fi^. The predictions from ResBos are shown as the red histogram and 
from ResBos with small-x broadening as the black histogram [which is visible principally in (e)]. 
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FIG. 3: (color online) Ratio of the corrected distributions of (1/V) x (da/d<f>*) to ResBos for: (a) \y\ < 1, (b) 1 < \y\ < 2 
and (c) 1 2/ 1 > 2. Statistical and systematic uncertainties are combined in quadrature. In (a) and (b) a \ 2 for the comparison 
of the dielectron and dimuon data, xfee/w) ' * s calculated assuming uncorrelated uncertainties. The yellow band around the 
RESBosprediction represents the quadrature sum of uncertainty due to PDFs (evaluated using the CTEQ6.6 NLO error 
PDFs QJj]) and the uncertainty due to the QCD scale (evaluated by varying the factorization and renormalization scales 
simultaneously by a factor of two). Also shown are the changes to the ResBos predictions when 172 is set to 0.66 (dotted 
blue line) and when the small-x broadening option is enabled (solid black line). 
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by the DO detector at the Fermilab Tevatron, we have 
studied with unprecedented precision the pj, distribution 
of Z/-f* bosons in dielectron and dimuon final states. In 
bins of boson rapidity, the normalised cross section is 
measured as a function of the variable </>*. Predictions 
from ResBos are unable to describe the detailed shape 
of the corrected data, and a prediction that includes the 
effect of small- x broadening is strongly disfavored. 

Tables of corrected (1/er) x (da/d(j)*) distributions for 
each \y\ bin and channel are provided [201 ] . 
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Supplementary Material 

Tables ITH4VII show the values of the corrected 
(1/cr) x (da I d(j)*) distributions for each \y \ bin and chan- 
nel. The first uncertainty is statistical and the second is 
systematic. The data are corrected back to the particle 
level, corresponding to the kinematic cuts and definitions 
of particle level leptons as described in the text. The in- 
tegral of (1/cr) x (da/dcj)*) for cj>* in the range 0-4.749 
is normalised to unity separately for each \y\ bin and 
channel, with the exception of the dielectron channel for 
\y\ > 2 in which the integral is normalised in the c6* range 
0-2.049 . 



TABLE II: The corrected (1/cr) x (da/dfi*) distribution for 
the dielectron channel and \y\ < 1. 

bin range (1/°") X (da / dtf)^) 



1 





.000-0 


.010 


13.242 ± 0.065 ± 0.020 


2 


0. 


.010-0. 


.020 


12.006 ± 0.062 ± 0.012 


3 


0. 


.020-0 


.030 


10.429 ± 0.058 ± 0.012 


4 


0. 


.030-0 


.040 


8.756 


± 


0.053 


± 


0.008 


5 


0. 


.040-0 


.050 


7.183 


± 


0.048 


± 


0.007 


6 


0. 


.050-0 


.060 


5.911 


± 


0.043 


± 


0.005 


7 


0. 


.060-0 


.071 


4.762 


± 


0.039 


± 


0.005 


8 


0. 


.071-0 


.081 


4.070 


± 


0.035 


± 


0.003 


9 


0. 


.081-0 


.093 


3.387 


± 


0.031 


± 


0.005 


10 


0. 


.093-0 


.106 


2.806 


± 


0.026 


± 


0.004 


11 


0. 


.106-0 


.121 


2.279 


± 


0.022 


± 


0.003 


12 


0. 


.121-0 


.139 


1.830 


± 


0.018 


± 


0.003 


13 


0. 


.139-0 


.162 


1.414 


± 


0.014 


± 


0.003 


14 


0. 


.162-0 


.190 


1.084 


± 


0.011 


± 


0.002 


15 


0. 


.190-0 


.227 


0.750 


± 


0.008 


± 


0.002 


16 


0. 


.227-0 


.275 


0.513 


± 


0.006 


± 


0.001 


17 


0. 


.275-0 


.337 


0.333 


± 


0.004 


± 


0.001 


18 





.337-0 


.418 


0.197 


± 


0.003 


± 


0.000 


19 


0. 


.418-0 


.523 


0.115 


± 


0.002 


± 


0.000 



bin range (1/cr) x (da/d<f>^) (xlOOO) 

20 0.523-0.657 61.731 ± 1.238 ± 0.090 

21 0.657-0.827 32.115 ± 0.798 ± 0.080 

22 0.827-1.041 16.496 ± 0.509 ± 0.071 

23 1.041-1.309 7.960 ± 0.324 ± 0.171 

24 1.309-1.640 3.882 ± 0.203 ± 0.087 

25 1.640-2.049 2.006 ± 0.133 ± 0.056 

26 2.049-2.547 1.068 ± 0.090 ± 0.033 

27 2.547-3.151 0.702 ± 0.067 ± 0.028 

28 3.151-3.878 0.389 ± 0.045 ± 0.015 

29 3.878-4.749 0.284 ± 0.036 ± 0.013 



9 



TABLE III: The corrected (1/ct) x (da / dcj>*) distribution for 
the dielectron channel and 1 < \y\ < 2. 

bin range (1/°") x {da/dcf)*^) 

1 0.000-0.010 14.235 ± 0.104 ± 0.053 

2 0.010-0.020 12.782 ± 0.099 ± 0.046 

3 0.020-0.030 11.035 ± 0.092 ± 0.030 

4 0.030-0.040 9.023 ± 0.083 ± 0.019 

5 0.040-0.050 7.268 ± 0.074 ± 0.011 

6 0.050-0.060 5.911 ± 0.067 ± 0.007 

7 0.060-0.071 4.988 ± 0.061 ± 0.006 

8 0.071-0.081 4.029 ± 0.053 ± 0.007 

9 0.081-0.093 3.346 ± 0.047 ± 0.006 

10 0.093-0.106 2.696 ± 0.040 ± 0.008 

11 0.106-0.121 2.267 ± 0.034 ± 0.010 

12 0.121-0.139 1.780 ± 0.027 ± 0.010 

13 0.139-0.162 1.308 ± 0.021 ± 0.007 

14 0.162-0.190 1.015 ± 0.016 ± 0.005 

15 0.190-0.227 0.683 ± 0.012 ± 0.005 

16 0.227-0.275 0.444 ± 0.008 ± 0.003 

17 0.275-0.337 0.282 ± 0.006 ± 0.003 

18 0.337-0.418 0.157 ± 0.004 ± 0.002 

19 0.418-0.523 0.089 ± 0.003 ± 0.001 
bin range (1/a) x (da/d<f>*) (xlOOO) 

20 0.523-0.657 45.296 ± 1.581 ± 0.371 

21 0.657-0.827 22.931 ± 1.019 ± 0.113 

22 0.827-1.041 10.886 ± 0.636 ± 0.092 

23 1.041-1.309 4.909 ± 0.378 ± 0.114 

24 1.309-1.640 2.848 ± 0.262 ± 0.073 

25 1.640-2.049 1.330 ± 0.161 ± 0.045 

26 2.049-2.547 0.921 ± 0.120 ± 0.044 

27 2.547-3.151 0.363 ± 0.069 ± 0.019 

28 3.151-3.878 0.226 ± 0.051 ± 0.010 

29 3.878-4.749 0.214 ± 0.048 ± 0.013 
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TABLE IV: The corrected (1/cr) x (da/d<f>*) distribution for 
the dielectron channel and \y\ > 2. 



bin range (1/°") x (do"/dc/>*) 



1 


000-0 010 


15.625 ± 0.361 ± 0.031 


2 


0.010-0.020 


14.288 ± 0.344 ± 0.023 


3 


0.020-0.030 


12.130 ± 0.319 ± 0.030 


4 


0.030-0.040 


9.514 ± 


0.281 ± 


0.027 


5 


0.040-0.050 


7.572 ± 


0.250 ± 


0.009 


6 


0.050-0.060 


6.311 ± 


0.226 ± 


0.016 


7 


0.060-0.071 


5.052 ± 


0.202 ± 


0.013 


8 


071-0 081 


3.991 ± 


0.175 ± 


0.008 


9 


0.081-0.093 


3.206 ± 


0.152 ± 


0.008 


10 


0.093-0.106 


2.533 ± 


0.126 ± 


0.006 


11 


0.106-0.121 


1.796 ± 


0.099 ± 


0.006 


12 


0.121-0.139 


1.658 ± 


0.087 ± 


0.006 


13 


0.139-0.162 


1.223 ± 


0.067 ± 


0.002 


14 


0.162-0.190 


0.767 ± 


0.047 ± 


0.005 


15 


0.190-0.227 


0.605 ± 


0.037 ± 


0.003 


16 


0.227-0.275 


0.378 ± 


0.025 ± 


0.003 


17 


0.275-0.337 


0.195 ± 


0.016 ± 


0.002 


18 


0.337-0.418 


0.114 ± 


0.011 ± 


0.001 


19 


0.418-0.523 


0.050 ± 


0.006 ± 


0.001 


bin 


range 


(1/cr) x {do/d<f)* v ) 


(xlOOO) 


20 


0.523-0.657 


24.457 ± 3.781 ± 0.616 


21 


0.657-0.827 


4.145 ± 1.382 ± 0.122 


22 


0.827-1.041 


2.906 ± 1.099 ± 0.148 


23 


1.041-1.309 


0.503 ± 0.356 ± 0.024 


24 


1.309-1.640 


0.154 ± 0.154 ± 0.014 


25 


1.640-2.049 


0.101 ± 0.101 ± 0.009 
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TABLE V: The corrected (1/cr) x (da/dcj>*) distribution for 

the dimuon channel and \y\ < 1. 

bin range (1/°") x {da/dcf)*^) 

1 0.000-0.010 12.992 ± 0.058 ± 0.027 

2 0.010-0.020 11.958 ± 0.055 ± 0.022 

3 0.020-0.030 10.263 ± 0.051 ± 0.015 

4 0.030-0.040 8.620 ± 0.047 ± 0.013 

5 0.040-0.050 7.115 ± 0.043 ± 0.010 

6 0.050-0.060 5.863 ± 0.039 ± 0.007 

7 0.060-0.071 4.868 ± 0.035 ± 0.006 

8 0.071-0.081 4.084 ± 0.031 ± 0.006 

9 0.081-0.093 3.400 ± 0.027 ± 0.005 

10 0.093-0.106 2.831 ± 0.024 ± 0.004 

11 0.106-0.121 2.320 ± 0.020 ± 0.003 

12 0.121-0.139 1.850 ± 0.016 ± 0.003 

13 0.139-0.162 1.439 ± 0.013 ± 0.003 

14 0.162-0.190 1.061 ± 0.010 ± 0.002 

15 0.190-0.227 0.779 ± 0.007 ± 0.002 

16 0.227-0.275 0.526 ± 0.005 ± 0.002 

17 0.275-0.337 0.331 ± 0.004 ± 0.001 

18 0.337-0.418 0.207 ± 0.003 ± 0.001 

19 0.418-0.523 0.117 ± 0.002 ± 0.001 
"bin range (1/a) x (da / 'd</>*) (xlOOO) 

20 0.523-0.657 66.026 ± 1.160 ± 0.479 

21 0.657-0.827 34.188 ± 0.753 ± 0.304 

22 0.827-1.041 17.056 ± 0.480 ± 0.256 

23 1.041-1.309 8.753 ± 0.315 ± 0.254 

24 1.309-1.640 4.774 ± 0.214 ± 0.130 

25 1.640-2.049 2.489 ± 0.139 ± 0.073 

26 2.049-2.547 1.394 ± 0.096 ± 0.044 

27 2.547-3.151 0.802 ± 0.066 ± 0.025 

28 3.151-3.878 0.536 ± 0.051 ± 0.025 

29 3.878-4.749 0.307 ± 0.034 ± 0.012 
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TABLE VI: The corrected (1/a) x (da/dcj>*) distribution for 
the dimuon channel and 1 < \y\ < 2. 

bin range (1/°") x {do/dcff^) 



1 





000-0 


.010 


13.404 ± 0.105 ± 0.051 


2 


0. 


.010-0 


.020 


12.008 ± 0.100 ± 0.026 


3 


0. 


.020-0 


.030 


10.647 ± 0.094 ± 0.016 


4 


0. 


.030-0 


.040 


8.755 


± 


0.086 


± 


0.013 


5 


0. 


.040-0 


.050 


7.230 


± 


0.077 


± 


0.016 


6 





.050-0 


.060 


5.804 


± 


0.069 


± 


0.019 


7 





.060-0 


.071 


4.972 


± 


0.063 


± 


0.011 


8 





.071-0. 


.081 


4.045 


± 


0.056 


± 


0.008 


9 





.081-0 


.093 


3.441 


± 


0.050 


± 


0.008 


10 


0. 


.093-0. 


.106 


2.820 


± 


0.043 


± 


0.009 


11 





.106-0 


.121 


2.330 


± 


0.036 


± 


0.006 


12 


0. 


.121-0 


.139 


1.824 


± 


0.029 


± 


0.004 


13 


0. 


.139-0 


.162 


1.414 


± 


0.023 


± 


0.003 


14 





.162-0 


.190 


1.066 


± 


0.018 


± 


0.002 


15 


0. 


.190-0 


.227 


0.756 


± 


0.013 


± 


0.001 


16 


0. 


.227-0 


.275 


0.514 


± 


0.009 


± 


0.002 


17 


0. 


.275-0 


.337 


0.326 


± 


0.007 


± 


0.002 


18 





.337-0 


.418 


0.200 


± 


0.005 


± 


0.001 


19 


0. 


.418-0 


.523 


0.107 


± 


0.003 


± 


0.000 



bin range (1/a) x (da/dcj)^) (xlOOO) 

20 0.523-0.657 54.251 ± 1.874 ± 0.638 

21 0.657-0.827 25.906 ± 1.162 ± 0.259 

22 0.827-1.041 12.306 ± 0.730 ± 0.261 

23 1.041-1.309 5.197 ± 0.439 ± 0.174 

24 1.309-1.640 2.536 ± 0.279 ± 0.085 

25 1.640-2.049 1.263 ± 0.206 ± 0.086 

26 2.049-2.547 0.503 ± 0.113 ± 0.038 

27 2.547-3.151 0.350 ± 0.104 ± 0.022 

28 3.151-3.878 0.209 ± 0.063 ± 0.014 

29 3.878-4.749 0.069 ± 0.042 ± 0.012 



